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We show that the conventionally defined partonic cross section for the production of isolated 
prompt photons is not an infrared safe quantity. We work out the case of e + e~ -> 7 + X m 
detail, and we discuss implications for hadron reactions such as pp — > 7 + X. 
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1 Terminology 

1.1 Photon Isolation 

In e + e~ and in hadron-hadron reactions at collider energies, prompt photons are observed and their 
cross sections are measured only if the photons are relatively isolated in phase space. Isolation is 
required to reduce various hadronic backgrounds including those from the electromagnetic decay 
of mesons, e.g., tt° — ► 27. The essence of isolation is that a cone of half-angle 5 is drawn about the 
direction of the photon's momentum, and the isolated cross section is defined for photons accompanied 
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by less than a specified amount of hadronic energy in the cone, e.g., 
denotes the energy of the photon. Instead of 5, the Fermilab collider groups use the variable R = 
a/ (A0) 2 + (A^) 2 , where A77 and Acft denote differences of rapidity and azimuthal angle variables. 
Theoretical predictions will therefore depend upon the additional parameters e/j and 5 (or R). The 
isolated cross section is not an inclusive cross section, and, as we discuss below, the usual factorization 
theorems for inclusive cross sections do not apply. Isolation removes backgrounds, but it also reduces 
the signal. For example, it reduces the contribution from processes in which the photon emerges from 
the long-distance fragmentation of quarks and gluons, themselves produced in short-distance hard 
collisions. 
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1.2 Conventional Factorization 

Much of the predictive power of perturbative QCD derives from factorization theorems 0. Con- 
ventional factorization expresses a physical quantity as the convolution of a partonic term with a 
nonperturbative long-distance matrix element. It requires that the partonic term, calculated pertur- 
batively order- by-order in the the strong coupling strength a s , have no infrared singularities, and 
that the long-distance matrix element be universal. Applied to the case of e + e~ — > AX, conventional 
factorization states 

(?e+e-^A+x(Q) = CT e +e'->c+x{x, Q/fl) D c ^ A (z, fj,) + O (J^J ■ (1) 

The intermediate state c may be a quark, antiquark, gluon, or photon. The symbol (8> denotes a 
convolution. Variable z = pa/Pc- The fragmentation function D c ->a(z, (J,) represents long-distance, 
small momentum scale physics; its dependence on the fragmentation scale fi is governed by the 
Altarelli-Parisi evolution equations. In situations in which the factorization theorem applies, the 
partonic hard-part ("the partonic cross section") o" e + e -_ >c+ x(^, Q/l 1 ) is infrared safe, i.e., finite when 
all infrared regulators are removed. 

Conventional factorization holds for inclusive prompt photon production e + e _ — * jX, demon- 
strated through next-to- leading order 0. However, since the isolated cross section is not an inclusive 
quantity, factorization need not hold and indeed does noti. Nevertheless, almost all existing calcula- 
tions of the cross section for isolated photon production assume its validity^. Following the standard 
calculational procedures of perturbative quantum chromo dynamics (QCD), we show that the par- 
tonic hard-part for the isolated photon cross section is not infrared safe. The infrared sensitivity 
shows up first in the next-to-leading order quark-to-photon fragmentation contributions. We use the 
terminology "breakdown of conventional factorization" to describe this result. 



2 Photons in e + e Annihilation, e + e — > 7 + X 

Electron-positron reactions offer a relatively clean environment for the study of prompt photon 
production in hadronic final states B. Since there are no complications from initial state hadrons, 
e + e~ — > 7 + X is a good process in which to examine QCD predictions in the final state, and the 
data may be a good source of information on quark-to-photon and gluon-to-photon fragmentation 
functions i 3i. In turn, these fragmentation functions are needed for predictions of photon yields in 
hadronic collisions. Hard photons in e + e~ processes arise as QED bremsstrahlung from the initial 
beams, radiation that is directed along angles near # 7 = and tt, and as final state radiation from 
direct and fragmentation processes. The topic of interest in this paper is the final state radiation. It 
populates all angles, with an angular distribution having both transverse, l + cos 2 ^ 7 , and longitudinal 
components H. 

Cataloging the contributions to e + e~ — > 7 + X, we may first list the lowest order partonic 
process: e + e~ — ► q + q, followed by quark or antiquark fragmentation into a photon, q — ► 7X. The 
0(a em ) direct contribution is represented by the three-body final state process e + e~ — > q + q + 7. 
Separation of the lowest order fragmentation and the 0{a em ) direct contributions is not unique; the 
fragmentation scale \x dependence relates the two. At 0(a s ) the contributions of interest come from 
the three-body final state processes e + e~ — > q + q + g, followed by gluon fragmentation, g — > 7A; 
and from e + e~ — > q + q + g, followed by quark or antiquark fragmentation, q — ► jX. In Fig. 1(a) 
and 1(b), we illustrate the set of Feynman diagrams used in the computation of the contribution 
from e + e~ — * q + q + g (with the implied understanding that the final quark or antiquark fragments 
into a photon). The set of real gluon emission diagrams in Fig. 1(a) results in an infrared sensitive 
contribution, associated with the region of phase space in which the emitted gluon momentum becomes 
soft. Likewise, the set of virtual gluon loop diagrams in Fig. 1(b) also results in an infrared sensitive 
contribution. In the case of inclusive photon production, the infrared singularities cancel once the 
results from the two sets are combined. However, in the isolated photon case, the restriction on the 
phase space accessible to the final state gluon in the set of Fig. 1(a) leads to an incomplete cancellation 




Figure 1: Feynman diagrams are shown in (a) and (b) for the 0(a s ) cross section ^J+ e -_ >qX ', both real gluon emission 
diagrams (a) (e + e~ — > qqg) and virtual gluon exchange terms (b) are drawn. In (c), an illustration is presented of an 
isolation cone containing a parton c that fragments into a 7 plus hadronic energy -Ef rag . In addition, the cone includes 

a gluon that fragments giving hadronic energy i? par ton- 

between the two sets. The 0(a s ) partonic hard part <r^ e - x (x,Q/ fi) is therefore not finite when 
the infrared regulator is removed. 

3 Breakdown of Conventional Factorization 

In this section, we summarize our explicit calculation that demonstrates infrared sensitivity of the 
partonic hard part for isolated prompt photon production in e + e~ — > jX. It is useful to examine 
the sketch in Fig. 1(c). In that sketch, we illustrate an isolation cone centered on a 7 produced 
through fragmentation of a parton c. For the fragmentation contributions, there are two sources of 
hadronic energy in the isolation cone: i) -Ef rag from fragmentation of parton c itself, and ii) -Epartons 
from final-state partons other than c that also happen to be in the cone. The total hadronic energy 
in the cone is -E£ad?ons = £f rag + partons- For an isolated 7, -E£ad?ons < E cut , where E cut denotes the 
arbitrary limitation on hadronic energy in the cone that is selected in experiments. We choose to write 
E C ut = ChE-y, an equation that defines the quantity e^. When the maximum hadronic energy allowed 
in the isolation cone is saturated by the fragmentation energy, E cut = £f rag , there is no allowance for 
energy in the cone from other final-state partons. In particular, if there is a gluon in the final state, 
the phase space accessible to this gluon is restricted. 

We make frequent use of the variable x 7 = 2Ey/y/s, and we define 

%crit = z ; • (2) 
1 + €h 

There are three cases Of lilt 6I*6St • X iy X CT'it 3 X?y X Cffif, 5 clIXCl X -y ^ > X Cfif, • Because of the isolation 

condition, the phase space constraints are different in the three regions. We summarize the physical 
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Figure 2: One-loop quark fragmentation contributions to the isolated and inclusive cross sections (a) as a function of 
x y — 2E 1 / y/s in e + e~ — » jX at y/s = 91 GeV, and (b) as a function of r in pp — ► 7X at ^/s = 1.8 TeV. 



situation in the separate regions and show that the next-to-leading order partonic term for quark 
fragmentation, Eqda™° e - x /d 3 pq, is infrared sensitive^ at and below the point x 7 = 1/(1 + eh)- 

3.1 Xl <1/(1 + e h ) 

When x-y < 1/(1 + e^), subprocesses with two-body final states do not contribute. This statement 
follows from energy conservation. In a two-body final state, S 7 +£ , f rag = \fs/2, and E^^ ons = -^frag = 
(1 — x 7 )y / i/2. Isolation requires Madrons — e hE-y = ^h x jV^/2. Correspondingly, the two-body final 
state processes contribute only for x 7 > x C rit- I n particular, since two-body final states are absent, 
there is no contribution from one-loop virtual gluon exchange diagrams. There is a contribution from 
the three-body final state real gluon emission diagrams. At the partonic level, these yield a well- 
known infrared pole singularity having the form 1/(1 — x q ) as x q = x-y/z — > 1. Because the one-loop 
virtual gluon exchange diagrams do not contribute, the infrared pole singularity remains uncanceled 

e+e —>qX 

If conventional factorization were valid, the fragmentation contributions to the physical cross 
section would be expressed in the factorized form 



E. 



^e+° e -^ 7 X _v- f 1 dz do^e—MsX / _ ^7 \ D c ^(z,5) 



dH 

J max\ x-y, — - — z d?p, 



- E r e+ ;;^ cX u c = ^ c ^ 7V ^^ ; (3) 



x c = 2E c /^/s. The sum extends over c = q,q and g. Function D c ^(z,5) is the nonperturbative 
function that describes fragmentation of parton "c" into a photon; its evolution is governed by the 
Altarelli-Parisi equations. The lower limit of integration results from the isolation requirement with 
the assumption that all fragmentation energy is in the isolation cone H. After convolution with 
D q ^(z), the inverse power infrared sensitivity of cr l e s ° e _^ qX , at the partonic level, yields a logarithmic 
divergence in the physical cross section (?™+ e - x proportional to En(l/x 7 — (1 + e^)). As shown in 
Fig. 2(a), this means that the isolated cross section would become larger than the inclusive cross 
section in the vicinity of x 7 — > 1/(1 + eh), a result that is certainly not physical. This infrared 
sensitivity in &™° e -_^ x signals a clear breakdown of conventional perturbative factorization. 

Our claim that conventional perturbative factorization breaks down has been contested!. How- 
ever, a careful reading of that paper shows that, after repeating our derivation, the authors come to 
the same conclusions as we do. To quote from the concluding paragraphs, "the [infrared] logarithms 



.. become large in the neighborhood of x 7 ~ x cr u "one has to study whether these [infrared] 
logarithms can be factored since they reflect long-distance effects"; and "they destroy the relevance 
of the perturbative expansion at least in the neighborhood of x cr it' > . 

3.2 x 7 = l/(l + e h ) 

When x 7 = 1/(1 + e/J, it is possible to have x q = x^/z = 1. The one-loop virtual gluon exchange 
diagrams contribute fully in this case, with contributions proportional to 5(1 — x q ). However, isolation 
constraints limit the phase space accessible to gluon emission in the real subprocess, e + e~ — > qqg. 
Consequently, the infrared divergences of the real and virtual contributions do not cancel completely 
in the isolated case, unlike the inclusive case. Working in n = 4 — 2e dimensions, after adding the 
real and virtual contributions, we find! 

E q | ^ + \ f | " Zn- \ \s (1 - x g ) + finite terms . (4) 

The presence of the uncanceled 1/e and 1/e 2 terms means that the regulator e cannot be set to 0. 
Therefore, at x q = 1, corresponding to x 7 = 1/(1 + e/J, the partonic term for quark fragmentation 
is infrared divergent, and the perturbative calculation is not well-defined. Conventional perturbative 
factorization again breaks down. 

In Fig. 2(a) we present the perturbatively computed one-loop quark fragmentation contributions 
to the physical cross sections for inclusive and isolated prompt photon production e + e~ — > 7 + X at 
y/s = 91 GeV. 

3.3 Relevance for Experiment 

The logarithmic divergence in the physical cross section da lso jdx^ for x 7 < x cr u is an integrable 
singularity. Data are presented in bins of finite width, related to the experimental resolution. If the 
perturbative divergence spans a very narrow region in x 7 , one that is smaller than the typical bin 
width, then the theoretical cross section will still be useful. For the situation examined in this paper, 
the perturbatively calculated isolated cross section exceeds the inclusive cross section over a region in 
x 7 that is not narrow. As indicated in Fig. 2(a), the region extends over an interval of about 4 GeV 
in Ey, larger than typical bin widths at LEP and SLC. This means that the perturbatively computed 
da tso /dx 1 cannot be accepted at face value. Conventional perturbative QCD leaves strong infrared 
sensitivity in the partonic hard-part at next-to-leading order, and we must look beyond fixed-order 
perturbation theory in order to derive an expression that makes physical sense. 

3.4 Origin and Possible Cure of the infrared Sensitivity 

For x 7 < x cr it, the infrared sensitivity of the form jj- ln( 1— ) comes from soft real gluon emission. 
An analogous case is the transverse momentum (qr) distribution for the production of the intermediate 
vector bosons W and Z or of massive lepton-pairs (the Drell-Yan process). At 0(a s ), the subprocess 
q + q —* 7* + g provides a singular distribution of the form a s /q T . Likewise, the thrust distribution 
in e + e~ annihilation, da/dT, is singular in the limit T — ► 1. Fixed-order calculations do not work as 
qx — ► or as T — ► 1, and resummation is invoked I 1 ™ 1 H. For isolated photon production, as discussed 
in this paper, a similar problem is encountered as x q — > 1, or, equivalently, as x 7 — ► • All- 
orders resummation of soft gluon radiation is therefore also suggested for isolated photon production. 
Owing to the limited phase space available for gluon radiation, a Sudakov suppression of the infrared 
divergence might be expected. 

While there are similarities, there are also significant differences between the isolated photon case 
and the qx distribution and thrust examples cited above. In the isolated photon case, the point of 
infrared divergence occurs within the physical region at a location determined by the experimenters' 
choice of e^, not at a fixed point at the edge of phase space (i.e., qr = 0, and T = 1). More 



importantly, in the Drell-Yan and thrust examples, the infrared divergence cancels exactly between 
the real emission and gluon loop diagrams. Absent in these cases are uncanceled 1/e 2 and 1/e poles 
that arise in the isolated photon case at x 7 = x cr u from the restricted phase space for real gluon 
emission, c.f., Eq. (||). In the Drell-Yan and thrust cases, the divergence is regulated, and soft-gluon 
resummation can be done in Fourier transformed impact parameter space. The uncanceled poles 
of the isolated photon case would be tantamount to an extra unbalanced 5(q T ) piece in the Drell- 
Yan case. For isolated photon production, resummation must be done at the partonic level, in the 
x q distribution, before the convolution is performed with the q — > 7A fragmentation function. A 
resummation procedure must be devised to handle the unbalanced 6(1 — x q ) problem of Eq. @. 

The presence of infrared sensitivity is a tip-off that non-perturbative effects are present and must 
be addressed. For example, in the example of the qr distribution in Drell-Yan case, non-perturbative 
functions gi are introduced in the implementation of resummation in impact parameter space 113. 
These unknown functions determine the behavior of the differential cross section at modest values 
of qT- The critical point x cr u of the isolated photon case, defined in Eq. is arbitrary since e^ is 
chosen experimentally. If e/j is very small, x cr a — > 1; there will be only a very narrow region over 
which resummation will matter, and fixed-order perturbation theory will be adequate. On the other 
hand, if we are interested in extracting fragmentation functions from the data, e^ must not be too 
small. There will then be a relatively large region in which non-perturbative functions analogous to g\ 
will play a significant role in fits to data. Their presence is a source of uncertainty for the extraction 
of quark-to-photon fragmentation functions in next-to-leading order. After resummation, instead of 
the divergence apparent in the solid curve in Fig. 2(a), the predicted isolated cross section below 
x 7 = x cr u will remain bounded from above by the inclusive prediction, and it will join smoothly at 
%crit to the form it takes in Fig. 2(a) above x cr a. 

3.5 Recapitulation for e + e~ — > 7 + X 

To recapitulate, in e + e~ — > 7 + X, the next-to-leading order partonic term associated with the quark 
fragmentation contribution is infrared sensitive when x 7 < 1/(1 + e/j). Conventional perturbative 
factorization of the cross section for isolated photon production in e + e~ annihilation breaks down 
in the neighborhood of x 7 = 1/(1 + e/j). The isolated cross section, as usually defined, is not an 
infrared safe observable and cannot be calculated reliably in conventional fixed-order perturbative 
QCD at, and in the immediate region below, x 7 = 1/(1 + e^). All-orders resummation offers a 
possible theoretical resolution, but this situation differs in important respects from other examples of 
successful application of resummation techniques. 

4 Hadron Collider Experiments 

In hadron collisions, A + B —* 7 A, we are interested in the production of isolated prompt photons 
as a function of the photon's transverse momentum, pj<. At next-to-leading order in QCD, one must 
include fragmentation at next-to-leading order. At this order, difficulties analogous to those in e + e~ 
annihilation are encountered also in the hadronic case. To illustrate the problem 0, we consider the 
contribution from a quark-antiquark subprocess in which the flavors of the initial and final quarks 
differ: q' + q' — * q + q + g, where q fragments to a 7. We keep only final state gluon radiation so 
that the results of our investigations of e + e~ annihilation can be exploited directly. We specialize to 
rapidity y 7 = and take equal values for the incident parton momentum fractions, x a = Xb = x = y/r. 
In the translation to the hadronic case, the variable x 7 becomes xt where xy = 2pt/v§ ~ xt/x with 
xt = 2pt/\/~s- The critical point at which infrared divergence is manifest in the e + e~ case tends to 
occur at relatively large values of x 7 . Owing to the correspondence x 7 ~ xt/x, the critical point 
in the hadron case occurs at small values of the parton momentum fraction x where its effects are 
enhanced by the parton densities. 

The special one-loop quark fragmentation contribution to the observed cross section takes the 



form 

-^7 — AB i-xT' X ~ f d/r $ 'o'(r) E~ — 9 iT 7 ^ (t) + other subprocesses . (5) 
d a l J x ^ HH d A i 

In Fig. 2(b), we show the integrand in Eq. (||) obtained after convolution with the parton flux <3?(t). 
We compare the integrands for the isolated and the inclusive cases, and we observe again that infrared 
sensitivity at fixed-order leads to the unphysical result that the isolated integrand exceeds the inclusive 
integrand. It is evident that the convolution with the parton flux substantially enhances the influence 
of the region of infrared sensitivity. 

The integrand of Eq. (||) is not a physical observable. Instead, the contribution to the hadronic 
cross section is the area under the curve in Fig. 2(b) from x\ to 1. The divergences above and below 
the point Xt = 1/(1 + (-h) [ or \fr = x t(1 + £h)\ are integrable logarithmic divergences, and thus they 
yield a finite contribution if an integral is done over all r. Indeed, after integration, the calculated 
isolated cross section may well turn out less than its inclusive counterpart since the unwarranted 
extra positive contribution associated with the logarithmic divergence in the region of small r may 
be more than compensated by the area between the inclusive and the isolated curves at large r. The 
issue is one of reliability. We stress that the perturbatively calculated one-loop partonic cross section 
E q da q ^,^ qX /d 3 p q , has an inverse-power divergence as x q — > 1 and has uncanceled 1/e 2 and 1/e poles 
in dimensional regularization i. This pole divergence for xt < 1/(1 + e/i) becomes a logarithmic 
divergence in E^da q 'g'^x /d 3 £ after the convolution with a long-distance q — ► 7 A fragmentation 
function. Although the logarithmic divergence near yfr = xt(1 + £h) is integrable, the isolated 
integrand should never exceed the inclusive integrand. It is not correct to accept at face value a 
prediction for an isolated cross section obtained from a perturbatively calculated integrand whose 
value exceeds that appropriate for the inclusive case (even if, after integration over r, the resulting 
isolated cross section is smaller than the inclusive). 

In the calculated isolated cross section, how large is the uncertainty associated with the infrared 
sensitivity of the next-to- leading order quark fragmentation terms? Referring to Fig. 2(b), we define 
the overestimate of the isolated cross section to be the area under the solid curve but above the dashed 
curve. We define the maximum isolated cross section to be the area under the lower of the solid and 
dashed curves. Doing so, we find that the next-to-leading order quark-fragmentation contribution 
to the physical isolated cross section is overestimated by 54%. While this overestimate is certainly 
large, one should bear in mind that the tree-level fragmentation contribution is not affected by the 
infrared uncertainty and that the tree-level term is larger than the one-loop quark fragmentation 
term that is of concern. Thus, the uncertainty in the overall fragmentation contribution will be 
typically only some fraction of 54%. Second, even though fragmentation may account for half or 
more of the inclusive prompt photon yield at relatively small values of pr, the fragmentation fraction 
is substantially reduced by isolation. We suggest therefore, that the net uncertainty in the physical 
isolated cross section associated with the next-to-leading order fragmentation terms may be at the 
10% level. 

5 Summary 

The results in both the e + e~ and hadronic cases challenge us to find a modified factorization scheme 
and/or to devise more appropriate infrared safe observables. 

Acknowledgments 

Work at Argonne National Laboratory is supported in part by the U. S. Department of Energy, 
Division of High Energy Physics, Contract No. W-31-109-ENG-38. 



References 



1. J. C. Collins, D. E. Soper, and G. Sterman, in Perturbative Quantum Chromodynamics, ed. 
A. H. Mueller (World Scientific, Singapore, 1989). 

2. E. L. Berger, X. F. Guo, and J. W. Qiu, Phys. Rev. D 53, 1124 (1996). 

3. E. L. Berger, X. F. Guo, and J. W. Qiu, Phys. Rev. Lett. 76, 2234 (1996); Phys. Rev. D 54, 
5470 (1996). 

4. P. Aurenche et al, Nucl. Phys. B 399, 34 (1993) and references therein; H. Baer, J. Ohnemus, 
and J. F. Owens, Phys. Rev. D 42, 61 (1990); E. L. Berger and J. W. Qiu, Phys. Lett. B 248, 
371 (1990) and Phys. Rev. D 44, 2002 (1991); L. Gordon and W. Vogelsang, P hys. Rev. D 
50, 1901 (1994); L. Gordon, Argonne report ANL-HEP-PR-96-60 Qhep-ph/961139lD , November, 
1996. 

5. ALEPH Collaboration, D. Buskulic et al, Z. Phys. C 57, 17 (1993) and references therein; 
DELPHI Collaboration, P. Abreu et al, Z. Phys. C 69, 1 (1995) and references therein; L3 
Collaboration, M. Acciarri et al, Phys. Lett. B 388, 409 (1996) and references therein; OPAL 
Collaboration, G. Alexander et al, Z. Phys. C 71, 1 (1996) and references therein. 

6. K. Roller, T. F. Walsh, and P. M. Zerwas, Z. Phys. C 2, 197 (1979); E. Laermann, T. F. Walsh, 
I. Schmitt, and P. M. Zerwas, Nucl. Phys. B 207, 205 (1982); E. W. N. Glover and A. G. Mor- 
gan, Z. Phys. C 62, 311 (1994); G. Kramer and B. Lampe, Z. Phys. C 34, 497 (1987) and Phys. 
Lett. B 269, 401 (1991); G. Kramer and H. Spiesberger, Proceedings of Annecy Workshop on 
Photon Radiation from Quarks, Annecy, France, December 1991, CERN report CERN 92-04, 
ed. S. Cartwright, pp. 26-40; P. Mattig, H. Spiesberger, and W. Zeuner, Z. Phys. C 60, 613 
(1993). 

7. E. L. Berger, X. F. Guo, and J. W. Qiu, in The Fermilab Meeting, Proceedings of the 7th 
Meeting of the Division of Particles and Fields of the American Physical Society, Fermilab, 
November, 1992; ed. C. H. Albright, P. H. Kasper, R. Raja, and J. Yoh (World Scientific, 
Singapore, 1993), pp. 957-960. 

8. ALEPH Collaboration, D. Buskulic et al, Z. Phys. C 69, 365 (1996); OPAL Collaboration, K. 
Ackerstaff et al, CERN report CERN-PPE/97-086 (July, 1997), submitted to Z. Phys. C. 

9. P. Aurenche, M. Fontannaz, J. Ph. Guillet, A. Kotikov, and E. Pilon, Phys. Rev. D 55, 1124 
(1997). 

10. J. C. Collins, D. E. Soper, and G. Sterman, Nucl. Phys. B 250, 199 (1985). For a recent 
treatment and list of references, see R. K. Ellis, D. A. Ross, and S. Veseli, Fermilab report 
FERMILAB-PUB-97/082-T (|hep-ph/970423g| ), April, 1997. 

11. M. Greco and Y. Srivastava, Phys. Rev. D 23, 2791 (1981); S. Catani, G. Turnock, B. R. Web- 
ber, and L. Trentadue, Phys. Lett. B 263, 491 (1991); G. Sterman, in QCD and Beyond, 
Theoretical Advanced Study Institute in Elementary Particle Physics (TASI'95), June, 1995, 
ed. D. E. Soper (World Scientific, Singapore, 1995). 



